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Abstract. The exact solution of the Dirac equation with the exact boundary condition 
is further investigated in the present paper. We introduce the unitary principle that can 
be used to disclose all logic paradoxes hidden in natural science and mathematics. By the 
unitary principle, we show that the existence of solution of differential equation requires 
an important modification to the constant of angular momentum in the original Dirac 
equation. Applying it to the hydrogen-like atoms, we derive the real solution involving 
the zero energy and the quantum proton radius as well as the neutron ionization energy. 
However, the expression for the discrete energy levels in the real solution only depends on 
the radial quantum number. How to treat the difference between the real expression for 
the discrete energy levels of the Dirac-Coulomb equation and the experimental observation 
of the atomic spectrums of hydrogen is not resolved yet. 



1. Introduction 

Darwin [1] and Gordon [2] first obtained the standard solution of the Dirac equation with a 
Coulomb potential in 1928. From then on, many different treatments to the Dirac equation 
were given one after another. In some literatures [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14], it 
constructs the second-order Dirac equations or the imitated first-order Dirac equation for 
the new upper and lower components by introducing the function transformation to the 
original Dirac equations. In recent years, we analyzed many different procedures to treat 
the Dirac equation and found that there were some substantial difficulties of mathematics 
hidden in those formal solutions[15, 16, 17, 18, 19, 20, 21]. For examples, the existence and 
uniqueness of solution for the bound-state case were destroyed there[22, 23, 24, 25, 26, 27]. 
Even if for the standard solution to the Dirac-Coulomb equation, one can also verify the full 
radial wave function of the S state at the origin not meeting the given boundary condition 
yet. In fact, the divergence of the wave function is one of the crucial logical difficulties for 
the relativistic wave equations. It should seek another treatment to the Dirac equation in 
order to obtain the real exact solution that satisfies all mathematical axioms and physical 
background. 

In 2006, we introduced the exact boundary condition with considering the radius of atomic 
nucleus instead of the traditional rough boundary condition without considering the radius 
of atomic nucleus and found another exact treatment [28] to the Dirac-Coulomb equation. A 
different exact solution, involving the different wave functions and the different expressions 
for the discrete energy levels, was hence obtained. However, by further analyzing the details 
of mathematical procedures of the various solutions under the traditional rough boundary 
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condition and our new solution under the exact boundary condition respectively, we found 
there was still another mathematical difficulty hidden in the original Dirac equation to 
distinguish the true from the false for its various formal exact solutions. It deals with why 
and how to construct the procedure to treat the constant of angular momentum. 25 years 
ago, in order to directly determine the ideas of right and wrong for all formal logic paradoxes 
hidden in theoretical physics, we summarized a general principle called the unitary principle. 
It provides great help to find the real solution to the Dirac equation. All of these works 
led to this paper. Here we try to come to hard mathematics logic to solve all mathematical 
difficulties hidden in various formal solutions of the Dirac equation, with the result that the 
exact quantum proton radius and neutron ionization energy were obtained. 



2. Traditional solution under rough boundary condition 

Here, it should point out a fact on quantum mechanics that once a wave equation with the 
condition for determining solution is given the rest would be pure mathematics only [29] . The 
preciseness of mathematics requires that one should not ignore any logical difficult hidden in 
a famous theory. If the bound state problems were put forward as the corresponding problem 
for determining the solution of equations of mathematical physics, one would always, first 
notice those fundamental differences of the formal solutions and the real solutions, thus 
the Dirac equation would have a different tendency to develop. For the Dirac equation 
describing the hydrogen-like atoms, the problem was finally translated into solving a system 
of radial differential equations for the upper and lower components of the original radial 
wave function 

1 fF(r)\ 
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obviously, the problem under the rough boundary condition is a pure mathematics problem 
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where c is the velocity of light in the vacuum, mis the rest mass of electron, a = e 2 / (AftEQhc) 
is the fine-structure constant, h = h/2n with the plank constant /land k = ±1, ±2, ±3, • • • 
is the artificial constant of motion. The conventional viewpoints have been that the above 
problem for determining solution only has the following solution [30, 31, 32, 33] 
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where n r = 0, 1, 2, • • • and the coefficient b v and d v satisfy the following system of recurrence 
relations 

mc 2 - E 

(s + v — k) b v + Zad v — ab v -\ d v -\ = 

{2A) m c* + E 

Zab v — (s + v + k) d v - b v -\ + ad w -i = 

ac 

where s = \J k 2 — Z 2 a 2 and a = ^/m 2 ^ — E 2 j he. 

(2.3) with (2.4) is often called the standard solution of Dirac-Coulomb equation., because 
the formula of the energy levels in (2.3) agrees with the experimental observation to the 
atomic spectrum of hydrogen, this standard solution has been regarded as succeeding in de- 
scribing the fin-structure of the hydrogen-like atoms. However, as k = ±1, R (r — > 0) = ±oo, 
the wave function for the S state at the origin does not meet the boundary condition, im- 
plying that the wave function in (2.3) should be a formal key to the problem for determining 
the solution (2.2). It usually explains the above mathematical paradox as that the Dirac 
wave function for the hydrogen ground state has a mild divergence at the origin. We con- 
sider that the "mild divergence" is not a law of nature! On the other hand, there exist other 
contradictions hidden in the standard solution of the Dirac-Coulomb equation such as the 
domain difficulty of the constant of motion. From 1988, we have tried to find an exacter 
treatment to the Dirac equation. It was found that the rough boundary condition in (2.2), 
which was first given in the Schrodinger equation, was not applicable to the Dirac equation. 



3. A DIFFERENT SOLUTION UNDER EXACT BOUNDARY CONDITION 

Usually, a theory does not accept any deduction concealing some mathematics errors, 
even though it agrees with the experimental observation. If the Dirac equation were the 
only exact equation for describing the hydrogen-like atoms, a correct method to derive the 
Dirac formula for the energy levels would exist. We should seek the correct treatment 
but not various specious explanations. In the different formal theories, in fact, the cited 
Dirac formulas for the energy levels are not the real deductions of the corresponding logic. 
However, the topic is often transferred to that the formula for the energy levels agrees with 
the energy spectrum of the Dirac equation. Such actions have become a huge resistance 
force so it is difficult to expose and criticize the errors of the formal logic in theoretical 
physics, even if it is completely wrong. Experiment is a main norm for verifying a physical 
deduction, but what if the experiment data did not deny some deductions derived from the 
incorrect theories? 

Malenfant made the contributions [34, 35] to overcome the divergence of relativistic wave 
functions at the origin. For the divergence of standard wave function of the Dirac equation 
for theS* state at the origin, it usually considers the fact that the atomic nucleus of the 
hydrogen-like atoms has finite radius. As one of the reliable treatments, it is generally 
assumed that equation (2.2) is correct only for r > 8 (= radius of proton). Inside 5, 
the potential has to be modified from a Coulomb 1/r potential to one corresponding to 
a spread-out charge distribution. So, to do a completely correct calculation, one solves 
the Dirac equation separately outside of and inside of 5, with two different potentials, and 
then matches, at r = S, the outside solution (i.e., the standard Dirac-Coulomb one) to the 
inside solution. The inside solution is finite at r = 0; its effect is to modify the energy-level 
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formula by a small correction that takes into account the finite radius of the proton. This 
idea just supports introducing the exact boundary condition to the wave equation. For 
the hydrogen-like atoms, by the exact boundary condition, it should change the standard 
problem for determining the solution (2.2) into the exacter problem for determining the 
solution to the Dirac-Coulomb equation outside 5 
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where 5 > is defined as the quantum radius of the atomic nucleus, it is an intrinsic 
parameter of the original Dirac equation. Using the variable substitution r = £ + 5, the 
exact boundary condition becomes R (£ — > 0) ^ ±oo, — > oo) = 0. Just like finding 
the standard solution to the problem for determining the solution (2.2), it easily obtained 
a different formal solution to the Dirac-Coulomb equation as follows [28] 
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where n r = 0, 1, 2, • • • , r > S and the coefficient b u and d v are determined by the following 
new coupled recurrence relations 
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This new form of solution of the radial Dirac-Coulomb equation is different from the stan- 
dard solution (2.3). First, the full wave function R (r) determined by (3.2) has not any 
divergence for all states, because it only holds for r > 5 > 0. It has relation to the quantum 
proton radius 5. Therefore, 5 is an intrinsic constant of the radial Dirac equation. Second, 
there is one big difference between the new expression for the discrete energy levels (the 
first formula in (3.2)) and the Dirac formula (the first formula in (2.3)). The new expres- 
sion for the discrete energy levels cannot describe the fine structure of hydrogen-like atoms. 
However, in mathematics, the above new solution is a necessary deduction of the original 
Dirac-Coulomb equation. In the standard solution of the Diarc equation, in order to obtain 
the Dirac formula for the energy levels , some mathematics rules had to be destroyed. Fur- 
ther studies show that the solution (3.2) is just the real form of mathematical deduction of 
the Dirac-Coulomb equation with the exact boundary condition. Of course, there, we show 
that it needs to determine the lvalue and make an important modification to the constant 
of motion by the existence of solution of differential equation. 
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4. Unitary principle and Real solution under exact boundary condition 

On the same Dirac wave equation describing the hydrogen-like atoms, the exact solutions 
obtained differ in different line of reasoning because of using the rough and exact boundary 
condition respectively. Two different forms of solution (2.3) and (3.2) form a logic paradox. 
One with a very good agreement with the experimental observation conceals the mathe- 
matical contradictions; another without any mathematical contradiction remains a sizable 
gap between the logic deduction and the experimental observation. It has been difficult to 
treat such a logic paradox in theoretical physics. 

In order to reach a consensus on treating those recognized theories concealing the math- 
ematical contradictions, 20 years ago, we generalized an important principle that is called 
the unitary principle. Describing the law of nature can choose different metrolo- 
gies, there are definite transforms among different metrologies, but the law of 
nature does not change per se because of choosing different metrologies. As 
the different mathematical forms in the different metrologies for describing the 
same law of nature are transformed into the same metrologies, it must be the 
same as the form in the present metrologies, 1=1, the transformation is unitary. 
This principle can be used to disclose all logic paradoxes hidden in science theory as well 
as mathematics and find the correct treatments to the corresponding problems. By the 
unitary principle, we ever used rough boundary condition to solve the second-order Dirac 
equations for the original components F (r) and G (r)of the wave function; it still derives 
the formal formula of energy eigenvalue for the hydrogen-like atoms. However, the formal 
wave function of the second-order differential equations is not equivalent to the formal wave 
function of the original first-order differential equations. In this logic connection, the first 
order Dirac equation is chosen as one metrologies, and the second order Dirac equation for 
the original radial components of wave function is chosen as another metrologies. In fact, we 
have proven that this equivalent second-order Dirac equation, all other second-order Dirac 
equationsand all kinds of imitated first-order Dirac equationshas not any real solution under 
the rough boundary condition. However, in the corresponding formal solution, for obtain- 
ing the Dirac formula for the energy levels , the boundary condition for the original Dirac 
equation and the uniqueness theorem of solution for the equation and so on were used to 
be destroyed. This case still appeared in the initial value problem of the Maxwell equation, 
seeing a paper in Chinese, in which as one of real solutions of the Maxwell equation the real 
mode of plane transverse electromagnetic was given [36]. 

According to the unitary principle, if the standard solution (2.3) were correct, the new 
exact solution (3.2) would deduce 5 = 0. However, by further investigating the existence of 
solution to the system of recurrence relations (3.3), we find that the system of recurrence 
relations (3.3) does not hold for any artificial constant of motions = ±1, ±2, ±3, • • • . This 
leads to the important modification to the radial Dirac-Coulomb equation, but the corre- 
sponding wave function and the expression for the discrete energy levels are still (3.2). It 
deduces the quantum proton radius Jand the neutron ionization energy. 

We focus on the existence of solution to the systems of recurrence relations (2.4) and (3.3). 
Recently, in writing other papers, we proved that in history constructing the constant of 
angular momentum for the radial Dirac equation is unreliable. This can be viewed from 
a hidden paradox in recurrence relations (2.4) for the standard solution (2.3). As n r = 0, 
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(2.4) becomes a system of equations 
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The first two equations are linearly dependent, so are the last two equations. By the 

second and third equations, one has (k + Vk 2 - Z 2 a 2 ^j = a^{mc 2 + E )/(mc 2 - E ) > 0, 

therefore k > 0. Consequently, the Dirac formula for the energy levels in (2.3) only holds 
for all positive numbers k = 1,2,3, ••• , it is inappropriate for constructing the constant 
of motion k = dbl, ±2, ±3, • • • in order to obtain the radial Dirac equation in the original 
Dirac theory. The existence of solution needs the modification to the constant of angular 
momentum by using the undetermined parameter C to replace the artificial constants, C is 
another intrinsic parameters of the Dirac equation. Thus, the problem for determining the 
solution to the Dirac-Coulomb equation outside 5 takes the following correct form 
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where C is determined by the existence of solution of the differential equation with the 
exact boundary condition itself. The above problem for determining the solution has the 
solution as the same as expressions (3.2), which is rewritten as follows 
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where n T = 0, 1, 2, ■ ■ ■ , r > 5. The expression for the discrete energy levels and the wave 
function is as same as (3.2), and the energy formula has no relation to the undermined 
parameter C. Replacing k by C, b u by f u and d u by g u in the recurrence relations in (3.3), 
it obtains the system of recurrence relations for determining the coefficients of the series f v 
and g u 
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The intrinsic parameter 5 and Co, Ci, C2,--- , C„ r corresponding to n r = 0, 1, 2, ■ • • 
will be uniquely determined by this system of recurrence relations, but only unique value 
of minimum quantum proton radius 5 is permitted. The expression for the discrete energy 
levels and wave function (4.3) with the recurrence relations (4.4) constitutes all real solutions 
of the radial Dirac equation for hydrogen-like atoms. 

5. Quantum proton radius and neutron ionization energy 

Seeing about the nonlinear homogeneous [37] system of recurrence relations (4.4), in which 
each of quantities E, C, f u , g u may depend on n. A sequence of numbers fi, f'2, • • • , f n 
andgi, <?2, ••• , 9n satisfying the system of recurrence relations (4.4) is uniquely determined 
once one of the values of fo and go, the so-called initial value, and the intrinsic constants 5 
and C n , are prescribed. According to the unitary principle, the values of 5 and C n will play 
a decisive role to affirm that the standard solution under the rough boundary condition is 
not the real solution of the Dirac-Coulomb equation. We now discuss the special solution of 
the system of recurrence relations (4.4) for zero energy state. In the real solution (4.3), as 
n r = 0, b\ = 62 = • • • = and di = (I2 = • ■ • = 0. In order to obtain the untrivial solution, 
order bo 7^ and do 7^ 0. Inserting n r = into the recurrence relations (4.4) yields 
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where a>o = \AngC 4 — Eq I he. The last two relations are equivalent to each other. By the 

first and second relations, one easily obtains 5 = he (C(j — Z 2 a 2 ) (2ZaEo)~ 1 . Substituting 
the first relation into the third or fourth relation yields 
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Substituting the second formula into the third formula in (5.1) gives 



(5.3) Co = Zay 

this formula is combined with (5.2) to read 



Eo 

Zah 



(5.4) E = 0, C = Za, 5 = 

mc 

Thus, as it can be seen, the constant of the angular momentum in Dirac equations should 
be C = C n with n = 0, 1, 2, • • • but not the fabricated constant k = ±1, ±2, ±3, • • • , and 
the quantum proton radius is not zero, 5^0. According to the unitary principle, the rough 
boundary condition is unsuitable for the Dirac equation describing hydrogen-like atoms and 
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the so-called standard solution (2.3) with (2.4) should be one of the formal solutions of the 
Dirac-Coulomb equations. Otherwise, it would read C = k and 5 = 0. 

The special energy eigenvalue Eq = corresponds to the unique energy state of neutron 
and neutron-like AE = E^ — Eq. As Z = 1, it deduces the neutron ionization energy and 
the quantum proton radius as follows 



(5.5) AE = mc\ 5 = - - 

where E^ is the energy of hydrogen-like atoms corresponding to n r = oo in the formula 
(4.3). The neutron ionization energy implies that the neutron is so hard that it would be 
broken up only by a photon with the energy rac 2 , and combining an electron and a proton 
into a neutron would emit a photon with the energy mc 2 . The term S = e 2 (47reo"ic 2 ) 
(=2.8117940285 fm) has been regarded as the classical electron radius. However, the above 
analysis shows that it should be the quantum proton radius. This is about triplication of 
the recent values that are respectively reported to be 0.88015 fm[38], 0.895±0.018 fm[39] 
and 0.89014fm, which the value is actually calculated by Lamb shift[40]. Admittedly, the 
accurate calculation (5.5) is the necessary deduction of the Dirac equation. The difference 
between the calculation and the result of indirect measurement strongly suggests the direct 
measurement to the proton radius. 

6. Summary and expectation 

As one of the important criteria for testing logic of nature science and mathematics, 
the unitary principle is introduced in this paper. We show that the existence of solution 
for the Dirac-Coulomb equation demands the correction for undetermined parameter C n . 
Replacing the rough boundary condition by the exact boundary condition and replacing the 
artificial constant k = ±1, ±2, • • • corresponding to angular momentum by strictly deriving 
the intrinsic undermined constant C n , are the mainly ameliorations to the Dirac theory 
for the hydrogen-like atoms. Being different from the traditional standard solution, the 
real exact solution of the Dirac equation is obtained by solving the Dirac equation with 
the exact (real) boundary condition, without any mathematical paradox. It predicts the 
quantum proton radius and the neutron ionization energy. The traditional standard solution 
is obtained by solving the Dirac equation with the rough boundary condition and the fictive 
constant of motion. It conceals some paradoxes that destroy the fundamental mathematical 
rules but have been covered up by those specious exploitations or mathematical distortion 
technique. Just like that the divergence of the Dirac wave function of the S state at 
the origin was called the "mild divergence". As one of the unreal solutions to the Dirac- 
Coulomb equation, the Dirac formula for the energy levels agrees pretty accurately with 
the experimentally observed hydrogen spectral. It seems to be a brainteaser! However, 
as the real solution of the Dirac Coulomb equation, the new formula for the energy levels 
does not give the fine structure splitting. In this case, out of question, it must choose the 
real solution and then make further correction. Otherwise, the physic would continually 
use much more wrong logic or false science to coin the anticipant deduction. We conclude, 
therefore, that the Dirac equation succeeds in describing the character of neutron combined 
by the proton and the electron. We know that the exact boundary condition will influence 
some conclusions of quantum field theory and quantum electrodynamics. Why cannot the 
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real solution of the Dirac equation with the exact boundary condition describe the fine 
structure of hydrogen atom? Does the Dirac equation contain another stricter deduction 
for energy eigenvalue having not been found? By the unitary principle, if the real solution 
presents no other evidence that the real wave function predicts results that agree with any 
other experimental data, then it would inspire new possibilities to look for the exacter 
amelioration to the Dirac equation. 
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